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Abstract 

For a given point set 5 in a plane, we develop a distributed algorithm to compute 
the a— shape of S. a— shapes are well known geometric objects which generalize the 
idea of a convex hull, and provide a good definition for the shape of S. We assume that 
the distances between pairs of points which are closer than a certain distance r > 
are provided, and we show constructively that this information is sufficient to compute 
the alpha shapes for a range of parameters, where the range depends on r. 

Such distributed algorithms are very useful in domains such as sensor networks, 
where each point represents a sensing node, the location of which is not necessarily 
known. 

We also introduce a new geometric object called the Delaunay-Cech shape, which 
is geometrically more appropriate than an q— shape for some cases, and show that it 
is topologically equivalent to a— shapes. 

1 Introduction 

Many applications call for detecting and tracking the boundary of a dynamically changing 
space of interest [1] [3]. We would expect any algorithm performing the task to include the 
following important properties: 1) the boundary output is geometrically close to the actual 
boundary, and 2) the interior of the boundary is topologically faithful to the original space. 
It is often the case that we are only given random samples from the space. We may then 
reconstruct the space by first placing balls of a certain radius around these points, and then 
by taking the union of these balls. A good exposition on relationship amongst the sampling 
density, the geometry of the underlying space, and the radius of the balls may be found in [TO] . 



In this paper, we start with the assumption that the union of the balls described above 
is a good approximation to the space of interest. Note that in some cases, this is by design. 
For example, in the case of systematic failures in sensor networks [3] , the failure in the nodes 
is caused by a spatially propagating phenomenon, and our aim is to track its boundary. In 
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this case, we construct a space by taking the union of balls of radius rc/2 around each node, 
where is its radius of communication. The radius of communication is the distance within 
which two nodes can communicate with each other. 

The problem may also be viewed as one of computing the boundary of a set of points, 
provided with some geometric information. Given the binary information, about nodes 
pair-wise location within a certain distance, distributive algorithms exist to compute such a 
boundary [2]. These algorithms are, unfortunately relatively slow, on account of the need 
for a global structure to reach a decision on the membership of a node or of an edge to the 
complex boundary. If on the other hand, we are provided with all pair-wise distances of 
nodes within a neighborhood, the above decision may be locally made by constructing an 
associated a— shape. 

Given a set of points S* in a plane, the a— shape introduced in |5j gives a generalization of 
the convex hull of S, and an intuitive definition for the shape of points. More importantly, 
an a— shape is the boundary of an alpha complex, which has the same topology as that of 
the union of balls. This relation amongst a— shape, alpha complex and the union of balls is 
contingent on certain relations between their parameters. We discuss this in detail in Section 
12.11 Such topological guarantees cannot be provided by the boundary computed in |2j. 

The Delaunay triangulation gives sufficient information to compute the alpha complex, 
and hence its boundary, the a— shape. If we only require the a— shape, less information 
would be necessary. The work in [7] shows that a global Delaunay triangulation is not neces- 
sary, and the alpha shape can be computed using local Delaunay triangulations. Given the 
edge lengths of a geometric graph, [7] constructs the local Delaunay triangulation by first 
building local coordinates and then computing the Delaunay triangulation. Computing the 
local coordinates, is however not robust and requires a high density of nodes for accuracy. 
When given the edge length information, we show that even local Delaunay triangulation is 
not necessary. 

When there is a sufficient density of nodes, computing local coordinates is accurate (prob- 
abilistically), and distributed algorithms exist for computing modified versions of Delaunay 
triangulation [HIS]. In this case, we define a certain Delaunay- Cech triangulation, which 
contains an alpha complex, and which we show to be homotopy equivalent. For bound- 
ary tracking-based applications, the boundary of Delaunay-Cech triangulation will serve as 
a better geometric approximation to the boundary, while preserving the topological features. 

Our contributions in this article are: 

• Given the distances between pairs of nodes whenever they are closer than rc > 0, we 
develop an algorithm to compute the a— shape for a range of parameters, where this 
range depends on r^. 

• We introduce the Delaunay-Cech triangulation, defined in Section 12. 2[ and show that 
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it is homotopy equivalent to the alpha complex. 

The remainder of the paper is organized as follows, in Section [21 we provide some back- 
ground information, along with a formulation of the problem. We describe the distributed 
algorithm for computing an a— shape in Section [3l The Delaunay-Cech triangulation and the 
Delaunay-Cech shape are defined in Section [2^21 while the proof of its topological equivalence 
to the alpha complex is given in Section |H We conclude in Section |5] with some remarks. 

2 Preliminaries 

2.1 Alpha complex and a— shape 

Consider a set of nodes V C M^, and a parameter r. Let Vi be the voronoi cell associated 
with node Vi & V in the voronoi decomposition of V. Define an alpha cell (a— cell) of Vi as 
a{vi, r) = ViH B{vi, r) where B{vi,r/2) is the closed ball of radius r/2 around Vi. The alpha 
complex, Ar (we are assuming V is implied in this notation), is defined as the nerve complex 
of the alpha cells, i.e., (fo,fi, . . ■ ,Vk) spans a fc— simplex in if Oi^i'^i) ^- Since the 
alpha cells are convex, the nerve theorem [SlIII] implies that the alpha complex has the same 
homotopy type as the union of the alpha cells, which in turn is equal to the union of the 
balls B{vi,r/2). 

Given a set of nodes C Q, and a parameter r > 0, the alpha shape, dAr, is a 
1-dimensional complex which generalizes the convex hull of V. To simplify the notation, we 
use {vi,Vj) to denote an edge in a graph, a 1-simplex in a complex or the underlying line 
segment. A 1-simplex {vi,Vj) belongs to dAr if and only if a circle of radius r/2 passing 
through Vi and Vj does not contain any other node inside it. By "inside" a circle, we mean 
the interior of the ball to which this circle is a boundary. We say that such a circle satisfies 
the ^^a— condition" . dAj. also contains all the nodes {vj} such that a circle of radius r passing 
through Vj satisfies the a— condition. 

For a 2-dimensional simplicial complex K, we define the boundary of K to be the union 
of all the 1-simplices (along with their faces), where each is a face of at most one 2— simplex, 
and all 0— simplices which are not faces of any simplex in K. The alpha shape dAr is the 
boundary of the alpha complex Ar [B]- 

2.2 Delaunay-Cech Shape 

For a set of nodes V C M."^ and a parameter r > 0, define the geometric graph Gr = (V, E) 
to be the set of vertices (V^) and edges (E), where e = {vi,Vj) is in E if the distance 
between Vi and vj is less than or equal to r. Let C{V,r) denote the Cech complex with 
parameter r (the nerve complex of the set of balls {B{vi,r /2)}) and let DT{V) be the 

tThe alpha shape is generally defined for points in M.^ for any dimension k. 
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Delaunay triangulation of V. We define the Delaunay-Cech complex DCr with parameter r 
as DCr = DT(y) nC(y, r). We will show in Section HJ that DCr is homotopy equivalent to 
Ar. We call the boundary of DCr, denoted by dDCr the Delaunay-Cech shape. 



3 Computing the alpha shape of points in M? 

In order to compute dAr, we take each edge in Gr and check if it is in dAr. If an edge 
e = belongs to dAr, then the length of the line segment is less than or equal 

to r. Otherwise, the a— condition cannot be satisfied. Edge e hence also belongs to Gr, and 
consequently, checking for all the edges in Gr is sufficient to compute dAr. 

Given an edge e = {vi,Vj), there are two circles of radius r passing through Vi and Vj. 



Let us call these circles C and C (see Figure 1(a)). The a— condition is satisfied if and only 



if at most one of C and C contains node(s) inside. 

We consider all the nodes in A/^ HA/^- (neighbors common to both Vi and vj) , and perform 
a series of tests to verify their location inside C and C. It is simple to see that considering 
nodes only in A/i fl f/j is sufficient. The diameter of both C and C is r. If Vk lies in one of 
the circles, the distance between Vk and either of Vi and Vj is less than r, and hence, Vk is a 
neighbor to both. 

We now derive the following: 

1. A test to see if a node lies in both circles C and C. This immediately determines that 
e does not belong to dAr. 

2. A test to see if a node lies in exactly one of the circles C and C. 

3. Given that there exists at least one node in one of the circles, a test to see if a subsequent 
node lies in the other. This also, immediately determines that e does not belong to 
OA. 

Let the angle subtended by the chord ViVj on the bigger arc (of either circle, see Figure 



1(a) ) be 6, hence making the angle subtended on the smaller arc ir — 6. The angle which the 
chord subtends at the center, u, may easily be computed using the law of cosines, and 9 is 
equal to uj/2. 

Let Vk G Mi nMj, and let ZviVkVj = (j)k. Then, ii (j)k > n — 9 , Vk lies inside both circles, 
and we immediately know that e does not belong to dAr. \i cf) < 9, it lies inside neither 
circle. Let Vk be the first node satisfying 9 < (j) < tt — 9. Then Vk lies in one of the circle. 
Without loss of generality, we assume Vk lies in C. 

Let vi be any subsequent node satisfying 9 < (pi < 71 — 9. If f ; is not a neighbor of Vk, then 
vi lies in C , since any two nodes inside a circle of diameter r will be neighbors. If Vk and vi 
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are neighbors, we know the length ||(ffc'y/)||. Using the law of cosines, we compute the angle 
^VkViVi which we call (3. If (3 = ZvkViVj + ZviViVj, vi lies in C, and if /3 = \ZvkViVj — ZviViVj\, 
Vi lies in C. Figure 1(b) demonstrates this relationship between the angles. 



computing the a— shape 



At each edge e = (fj, Vj) in G, 
compute 9 

for each Vk E MiH Afj 
compute 0fc 
if (pk > 'n' — 

e ^ dA, terminate 
if 0fc < 0, 

continue to next node 
ife<(l)<n-e 

is Vk the first node satisfying this condition? 

assign Vk to C 
else 
compute (3 

if (3 = \ZvkViVj - ZviViVjl 

continue to next node 
else 

e ^ dA, terminate 
eedA 

Table 1: Algorithm for computing the a— shape. Note that all the computations require 
only local information. 



The algorithm terminates when we determine that both circles C and C contain at least 
one node, or there are no more nodes in A/^ fl N'j to consider. In the former case, the edge 
e does not belong to dAr and in the latter, e belongs to dAr- Clearly, we can use the 
same algorithm to compute the a— shape for any parameter < g < r. The algorithm is 
summarized in Table [H Figure [2] shows dAr for a set of points in computed using the 
algorithm in Table [T] The shaded region in the Figure is the union of balls of radius r/2 
centered at each point. Note that the a— shape is a boundary of an object which is homotopy 
equivalent to the shaded region. 



4 Relation between DCr and 

Consider the Delaunay-Cech complex DCr as defined in Section I2l2| and the a— complex Ar, 
as defined in Section [2?T1 We will show that DCr has the same homotopy type as Ar, by show- 
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(a) (b) 

Figure 1: (a) shows the cycles C and C which pass through vi and Vj. (pk satisfies 6 < 
(f) < 71 — 6 and hes in C. (b)shows the angle relationships. When vi also lies in C, then 

/3 = \ZvkViVj — ZviViVjl, and if vi lies in C, (3 = ZvkViVj + ZviViVj 




Figure 2: a— shape with parameter rc/2 for a set of points in computed using algorithm 
in Table [1] The shaded region is the union of balls of radius rc/2 centered at each point. 

ing that there exists a bijective pairing between the 1-simplices and 2-simplices in DCr \ Ar, 
such that the pairing describes a homotopy collapse. Note that both A^. and DCr do not 
contain any simplices of dimension greater than 2. Figure |3] shows such a homotopy collapse. 




Figure 3: Homotopy collapse of an edge into an adjacent 2-simplex 
Let F{Gr) be the fiag complex of Gr- Define the complex Ur as Ur = DT(y) fl F{Gr). 
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(a) r e r/=/' (b) r e Te"/' 

Figure 4: Examples of the triangle sets Tl"^'^ and 77^^^. The circles in (a) have are centered 
at the nodes with radius r/2. Since the circum-radius is less than r/2, the circles have a 
common intersection with positive area. 

Since the Cech complex CiV^r) is a subcomplex of F{Gr), DCr is a subcomplex of Ur- 

Let Te denote the set of all 2-simplices to which e is a face in Ur, TT"^^ ^ Te denote the 
2-simplices in % with circum-radius less than or equal to r/2, and iT^'^ ^ % denote the 
2-simplices in Te such that the angle opposite e is greater than tt/2. Figure H] shows examples 
of triangles with these properties. 

In order to show the existence of a paring, we first analyze the triangles surrounding an 
edge e G Ar). The following lemma characterizes the 2-simplices in Ar in terms of their 
circum-radius. 

Lemma 4.1 A simplex (f i, f2, v^) in Ur is in Ar iff the circum-radius of the triangle (f i, f2, ^3) 
is less than or equal to r/2. 

Proof Since Ur C DT{V), {vi,V2,V3) G Ur ^ ("^i, "^2, "^3) ^ DT(y). The circumradius is 
less than or equal to r/2, iff circumcenter belongs to all a— cells a{vi), a{v2) and 0(^3). 
This results in the three a— cells having a non-empty intersection, hence {vi,V2,vs) G Ar. 
■ 

Lemmas 14. 21 and 14. 31 together impose conditions on cardinality of the sets 77"^'^ and T^^'^. 
We utilize these conditions in Lemma 14.11 to show the existence of the pairing. 

Lemma 4.2 Denote by rUe the midpoint of the 1-simplex e E Ur- rrie is a witness for e iff 

t:" = 

Proof nie is a witness for e iff there does not exist any other node inside the circle with e 



as the diameter (illustrated in Figure 5(a)). This occurs if and only if the angle opposite e 
in any incident triangle is acute. ■ 
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(a) Construction for Lemma 



(b) Construction for Lemma 



Figure 5: 

Lemma 4.3 Consider the following statements 

• Si : 17'^ = 

• S-i'. e E Ar 

Si y S2 is a necessary and sufficient condition for S3. 

Proof Let e = (fi,f2)- For sufficiency: from Lemma \A.2\ Si implies rrie is a witness for Vi 
and V2, and S2 implies 3 a witness (the circumcenter of one of the triangles in 77^"^) ■ For 
necessity: if 5*3 is true, then there exists a witness for e. If me is a witness, then 5*1 is true. 
If me is not a witness, then e shares a witness with a 2-simplex which is in A^. From Lemma 
14. H this implies 7^ 0. Therefore, 5*2 is true. ■ 



The above Lemma suggests the existence or non-existence of types of triangles surround- 
ing an edge in Ar. Lemma and Theorem 14. II further refine this relationship, and precisely 
identify the triangle to be removed when an edge is removed from DCr- 

Lemma 4.4 7/77^^ ^ 0, then \17'^\ = 1. 

Proof Suppose iT^"^ 7^ 0. Let r G iT^'^, and let the angle opposite e in r be with > 7r/2 
(see Figure 5(b) ). Let ri 7^ r be incident on e, with v being the opposite vertex. Since r G Ur, 
V does not lie inside the circum-center of r. This implies that the angle opposite e in ti is 
less than vr — which is less than tt/2. ■ 

Let Ck{K) denote the fc— simplices in the complex K. 



^we use the notation (vi, W2, ^3) to denote both the simplex and the underlying triangle 
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(a) Ar^/2{V) (b) DCr^/2{V) (c) Ar^/2{V) super-imposed over 



Figure 6: Figure shows the homotopy equivalence between Ar^/2{V) and DCr^/2(y). The 
shaded region is Rc. Note that DCrc/2(y) is a better geometric approximation to Rc than 

Theorem 4.1 Let Tr = C2 {DCr) \ C2 {Ar) and 8r = Ci {DC,) \ Ci {Ar). There exists a 
bijective pairing P : £r ^ Tr such that e is a face of P{e). 

Proof Let e G DCr but e ^ A^, from Lemma 77"^'^ = and TT^"^ 7^ 0- Since we assume 
e ^ Ar, e cannot be a face of any 2-simplex in Ar. Owing to the condition 7^"^^ = 0, Lemma 
14.11 ensures that this is indeed the case. Also, from Lemma [4.4[ \TT^'^\ = 1. Let r G TT^'^- 
Note that r is unique, and r ^ Ar. Further, since r is an obtuse triangle, r G Ciy,r), and 
this implies r G DC*,.. The pairing P is then defined as P{e) = r. ■ 

For any simplicial complex K, let ai and (T2 be simplices of dimension 1 and 2 such that 
0"! is a face of o"2. Then, there exists a deformation retraction Fg-^ : i^' ^ \ (cti U Int{a2)), 
which "collapses" ai into (T2. Therefore, K is homotopy equivalent to \ (cxi U Int[a2)). 

The removal of edges and triangles Tr describes a finite sequence of deformation 
retractions via the pairing P. When we collapse all the edges into their paired triangles, 
the resulting complex is Ar. Each collapse is a homotopy equivalence, and a composition of 
homotopy equivalences is a homotopy equivalence. This leads us to our main theorem: 

Theorem 4.2 The complexes DCr O'nd Ar are homotopy equivalent. 

Figure [6] illustrates the above theorem using an example. Note that Ar and DCr are 
homotopy equivalent to each other and both are homotopy equivalent to Rc (the shaded 
region). Further, as seen, DCr is a better geometric approximation to Rc than Ar. This is 
simply because Ar is a sub-complex of DCr. 

5 Conclusion 

The algorithm described in Section [3] takes the edge lengths as inputs and outputs the alpha 
shapes. We make no further assumptions on the node density, and we need not compute any 
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coordinates. The decision about an edge belonging to an a— shape is carried out by only 
looking at the local information, i.e., considering only the points within a certain distance, 
and may therefore be implemented distributively. In Section \2.2\ we define the Delaunay- 
Cech complex which contains the alpha complex. Its boundary, defined as a Delaunay-Cech 
shape, is therefore a better geometric approximation for the union of balls with an appro- 
priate radius. We also show in Section H] that, like the a— shape, the Delaunay-Cech shape 
remains topologically faithful to the underlying space. 
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